The power spectrum of the atomic dipole moment 
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This paper calculates the power spectrum S(u) of the electric field generated by the atomic dipole 
moment of a laser-driven two-level system from an open quantum systems perspective. Its shape 
is similar to the shape of Mollow's resonance fluorescence spectrum but there are some differences. 
For sufficiently strong laser driving, there are two sidebands but their relative height is reduced. 
Moreover, the amplitude of this spectrum has a different dependence on the laser Rabi frequency 
Q. It does not vanish when Q tends to zero. The calculation of the spectrum which we present 
here involves less approximations than the calculation of Mollow's spectrum and constitutes an 
interesting alternative property. 
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I. INTRODUCTION 

Recent years have seen a large number of experiments 
which study light-matter interactions at the quantum 
level using trapped ions p~j , single quantum dots [2] , color 
centres [3], and molecules on surfaces [TJ. What all of 
these systems have in common is that they enable the 
excitation of a single trapped electron into a well-defined 
excited state which can result in the spontaneous emis- 
sion of a photon. So-called artificial atoms, like quantum 
dots and color centres, often have much stronger spon- 
taneous decay rates T and much smaller transition fre- 
quencies Wo, thereby providing a new testing ground for 
quantum optical standard approximations. 

Moreover, recent technological advances [5] opened 
new possibilities in measuring extremely small forces. For 
example, Usenko et al. [6] recently demonstrated the de- 
tection of subattonewton forces at millikelvin tempera- 
tures by using a superconducting quantum interference 
device. Other authors employ quantum point contacts 
as sensitive displacement detectors in high precision ex- 
periments which aim at quantum limited displacement 
detection |U [TUj . The technology needed to combine a 
single quantum dot and an atomic force microscope can- 
tilever in a single experimental setup has already become 
available [IT] . 

Motivated by these developments, this paper calculates 
the power spectrum S(ui) of the electric field generated by 
a laser-driven atomic two-level system. We theoretically 
analyse the signal f(t) which is given by the expectation 
value (E(x))t of the electric field generated at a time t 
at a detector position x. Due to constant environmen- 
tal resetting, the free radiation field contains on average 
less than a single photonic excitation [T2TfTo] . Since sin- 
gle photons have a zero electric field expectation value, 
the only contribution to the signal f(i) comes from the 
atomic dipole moment r. Since the atomic dipole mo- 
ment of an atomic two-level system is essentially given 
by the Pauli operator er x = a + + er~ , we find that 

/(*) = K)* (i) 

up to an overall factor which does not depend on the 



atomic dynamics. As we shall see below, the power spec- 
trum S(us) associated with the signal in Eq. (fTl) is a direct 
measure for the coherence of the atomic source. Like the 
Mollow spectrum, it exhibits three peaks with the central 
peak being at w = ujq. The main difference is a signifi- 
cant reduction of the relative height of the sidebands and 
a different dependence of the amplitude of S(ui) on the 
laser Rabi frequency Q. 

The power spectrum of the signal /(£) = (E(x)) t pro- 
vides a direct measure of the atomic dynamics, but its 
derivation does not involve as many approximations as 
are used in conventional derivations of Mollow's fluores- 
cence triplet [TfJ>lU9| . For example, we treat the atomic 
system and the surrounding free radiation field as an 
open quantum system and assume a continuous loss of 
memory between the atomic system and its environment 
[2TJ] . The first measurements of the Mollow triplet were 
performed in the seventies [2"fll24| . Recently there has 
been a huge revival in this type of experiment. Quan- 
tum systems under investigation include atoms [25] . sin- 
gle molecules [T], quantum dots [55] [27], and quantum 
dots inside an optical cavity [5S1 HH] . 

There are four sections in this paper. In Section [il] we 
discuss the theoretical background for our calculations. 
We introduce the electric field operator E(x), summa- 
rize the atomic time evolution of a laser-driven atomic 
two-level system, and shortly review Mollow's spectrum. 
In Section we derive the power spectrum of the atomic 
dipole moment and compare its features with the prop- 
erties of Mollow's spectrum. Finally, we summarise our 
findings in Section [IV[ 



II. THEORETICAL BACKGROUND 

Quantum optical systems with spontaneous photon 
emission usually need to be described as open quantum 
systems. For example, when analysing an atomic sys- 
tem coupled to a free radiation field, it is assumed that 
a photon-absorbing environment resets the free radiation 
field on a coarse grained time scale At onto its vacuum 
state such that, at any given time t, there is at most one 



photon in the kA-modes of the free radiation field [T2]- 
H5l [31] . The vacuum is the einselected state [30] of a free 
radiation field with a photon- absorbing environment. In 
the absence of a photon source, it is the only state which 
does not evolve in time. This makes it the preferred state 
into which the radiation field relaxes rapidly once a pho- 
ton has been placed into it. Based on this assumption, 
the standard Markovian master equation for atomic sys- 
tems can be derived without having to employ a wide 
range of approximations [31] . 

The most important of these approximations is the 
rotating wave approximation. It can be avoided by 
means of a unitary transformation producing a diago- 
nal Hamiltonian for which the bare photon vacuum with 
the source (e.g. the atomic system) in its ground state is 
the true ground state of the total Hamiltonian [35]. Only 
within this rotating wave representation, continuous en- 
vironmental resetting onto the photon vacuum yields a 
zero stationary state photon emission rate, regardless of 
whether or not the Markovian approximation of extend- 
ing time integrals to infinity has been made [31] . Since 
it is physically reasonable to assume that the stationary 
state photon emission rate of the source be zero in the 
absence of any external driving we identify the vacuum of 
the rotating wave representation as the preferred photon 
vacuum onto which the environment resets. By express- 
ing a physical observable in terms of the photon creation 
and annihilation operators of the rotating wave repre- 
sentation we can then determine the effect of continuous 
environmental reseting on its dynamics. 



A. The electric field operator 

The signal /(£) which we analyse in this paper is the 
expectation value of the electric field scattered by a laser- 
driven atomic system. In the rotating wave representa- 
tion the relevant observable E at a detector with position 
x is the sum of two components [31 , 
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where a kA is the annihilation operator (in the rotating 
wave representation) for a single photon with wave vec- 
tor k, polarisation A, and frequency uj k , while eo and V 
are the dielectric constant and the field quantisation vol- 
ume. Since single photons have vanishing electric field 
expectation values all non-zero contributions to (E(x)) f 
come from this second term in Eq. p]) provided we as- 
sume the environment continuously resets the radiation 
field into the rotating wave vacuum. 

Neglecting an overall factor which depends only on pa- 
rameters like the actual position x of the detector but 



not the state of the atomic system and taking into ac- 
count that the atomic dipole moment r is proportional 
to cr x = er + + o~ in the Schrodinger picture, we arrive 
at Eq. (fTl). We remark that if overall factors are to be 
ignored, then using an expression for the electric field in 
another gauge, e.g. the Coulomb gauge, also results in 
Eq. p]) provided it is assumed that the radiation field 
is reset onto the photon vacuum belonging to that rep- 
resentation. Thus, when overall factors are ignored, the 
only nonzero contribution to the electric field expecta- 
tion value is equal to the expectation value of the static 
longitudinal electric field. We emphasise however that 
if the environment continuously resets the radiation field 
into the rotating wave vacuum, then what would actually 
be measured is the component represented by the second 
term in Eq. (pj which is not merely the longitudinal elec- 
tric field. 



B. Atomic master equations 

To calculate (<7 X )*, we employ in the following the stan- 
dard quantum optical master equation of a resonantly 
driven atomic two- level system with ground state |1), ex- 
cited state 1 2) and spontaneous decay rate T, 
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Here a + = |2)(1| and o~ = |1)(2| are the atomic lowering 
and raising operator, Hujo denotes the energy difference 
between the atomic levels, and f2 is a real laser Rabi fre- 
quency. Moving into the interaction picture with respect 
to 



Hq = huiQ o + a 



(4) 



the master equation in Eq. ^ becomes time- independent 
and can be solved analytically. Doing so and using the 
notation p id = (i\p A i\j) with p 22 = 1-pn and p 21 = p* 12 , 
we obtain 



Re pM*) = e- r */ 4 Rep 12 (0), 
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and the stationary state density matrix elements of the 
atomic interaction picture density matrix p^i given by 
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and 



Repf 2 = 0. 



(8) 



(9) 



The atomic system under consideration possesses a sta- 
tionary state only in the interaction picture. In the 
Schrodinger picture, its off-diagonal matrix elements be- 
come time-dependent. As a result we choose to remain 
in the interaction picture wherein the observable er x , in 
which we are interested, becomes time-dependent. 



C. The Mollow spectrum 

As an alternative to the power spectrum considered 
here, laser-driven atomic two-level systems are often 
characterised by the so-called Mollow triplet of resonance 
fluorescence [TTJ [TB]. Since the spectrum we consider 
has many similarities with Mollow's spectrum, we briefly 
summarise the main characteristics of the latter. It is de- 
fined as the Fourier transform of the correlation function 
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(10) 



Calculating this expression with the help of the master 
equations presented in the previous section yields }18| 
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with n defined as in Eq. ftfh, the detuning S given by 

S = UJ — U!q , (12) 



and with 



nvioi 



2ft 2 



T 



ir 



4(r 2 + 2ft 2 ) 16/x 



12ft 2 



r 2 + 2ft 2 



(13) 



This spectrum has its maximum at the atomic transition 
frequency ujq. In the case of sufficiently strong driving, 
two sidebands of equal height appear at u = luq — \i and 
ui = ujq + /i. For weak driving, the sidebands vanish and 
there is only a single peak. 

Unfortunately, there seems to be no clear operational 
definition of the power spectrum defined by the correla- 
tion function in Eq. (10 1. In order to relate this quan- 



tity to the properties of photons emitted by the atomic 
system one starts with the assumption that the atom- 
field system can be viewed as closed. One then solves 



for the dynamics of the electric field within the Heisen- 
berg picture and employs the rotating wave approxima- 
tion to simplify the result. In this context the rotating 
wave approximation constitutes the identification of the 
positive frequency component of the electric field (pho- 
ton annihilation operator) with the atomic lowering op- 
erator. Similarly one identifies the negative frequency 
component (photon creation operator) with the atomic 
raising operator. Thus, within the rotating wave approx- 
imation the solutions of the equations of motion for the 
positive and negative frequency components of the elec- 
tric field are effectively of the form E T (t) ex a ± (t) , with 
the result that Glauber's first order correlation function 
(E~(£ + r)E + (i)) gives Eq. (10) when it is assumed that 



the atom is in the stationary state. The correlation func- 
tion (E~(£ + t)E + (£)) is supposed to furnish a measure of 
correlations in photon number at two distinct times |34j , 
but its exact physical interpretation is not entirely clear 
|35j . Nevertheless its Fourier transform is usually taken 
as the definition of the fluorescence spectrum. 



III. THE POWER SPECTRUM OF THE 
ELECTRIC FIELD 

In classical physics, the power spectrum S(ui) of a sig- 
nal /(£) is defined as the modulus squared of its Fourier 
transform f(ui) which implies 

5(«) = l/MI 2 . (14) 

For a real signal /(£), this equation yields 

/•OO 

dt / di'e- 1 ^- 4 ') /(*)/(*')• (15) 
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After substituting r = t — t' , one can easily see that the 
power spectrum of a signal is the Fourier transform of its 
two-time correlation function G(t), 

/oo 
dre- iwT G(r) (16) 

-OO 

with G(r) defined as 

G(t) = / dtf(t)f(t + r). (17) 



In quantum physics, the correlation function G(t) can 
reflect a very rich inner dynamics of the system under in- 
vestigation. A measurement at time t changes the state of 
a quantum system in a non-trivial way. As a result, mea- 
surement outcomes can be highly correlated with previ- 
ous measurements. 

Let us now have a look at the power spectrum of the 
atomic dipole moment of a laser-driven atomic two-level 
system which is given by /(£) in Eq. pj. If the signal /(£) 
reaches a stationary state / ss , then fjt) = / ss most of the 
time and G(t) in Eq. JT7] simplifies to G{t) = /(0) /(r) 
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FIG. 1: The normalised power spectrum S(lj)/S(ljq ) an d the 
normalised Mollow triplet S'moi(w)/S'moi(^o) in Eqs. ( 24 1— ( 13 1 
as a function of u for T = 10 8 Hz, cj = 10 15 Hz, and n = 4I\ 
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FIG. 2: The normalised power spectrum S(lu)/S(ujo) and the 
normalised Mollow triplet Sm i(oj) / Smoi(^>o) as a function of 
oj for the same F and ojq as in Fig. IT] but for Q = 0.5 F. 



if at t = the system is in its stationary state. How- 
ever, the signal we consider evolves constantly in time 
and therefore, we cannot calculate G(t) via this simplifi- 
cation. We can, however exploit the fact that the system 
possesses a stationary state and calculate G(r) using the 
interaction picture with respect to H in Eq. Q. Within 
this interaction picture, the observable a* becomes time- 
dependent and equals 



<r x (t) = a + e' 1Uot + 



Joj t 



(18) 



The time-independent eigenvectors of this operator, 

|A li2 (t)) = (|1)± 6^*12))/^ (19) 

oscillate constantly in time. Moreover, one can show that 
its eigenvalues are given by 



Ai. 2 = ±1 



(20) 



The equation implies that when detecting the electric 
field E(x), there are only two possible measurement out- 
comes. The laser-driven atomic system exerts a force of 
a certain size onto a test charge which points in one out 
of two opposite directions. 



Using the notation in Eqs. (19) and (20) and normal- 



ising the correlation function G[t) in Eq. (17) as usual, 
we find that G(t) equals 



bf T dt £ A ^ 
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with the superoperator T T summarizing the atomic time 
evolution in Eqs. ^ and (J6|; 



T T (pAl(t)) = pAl(t + T) , 



(22) 



and with JPj(t) being the projector onto the time- 
dependent state \Xi(t)), 



Pi(t) = \Xi(t))(Xi(t)\, 



(23) 



In other words, G(t) is the expectation value for the 
product of the results of two subsequent measurements 
of the atomic dipole moment <7 X - one taken at t and 
one taken at t + r. The state immediately after the first 
measurement is the eigenstate corresponding to the first 
measurement's outcome. It is independent of the initial 
state of the system (at t — 0) due to the averaging over 
an infinitely long time scale (T — > oo). This implies 
PAi(t) = Pai a * almost all times t. 



Taking this into account and evaluating Eq. (21) with 
the help of Eqs. (|Bj>— @, yields 
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with the coefficients j3± given by 
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Substituting this result into Eq. 



(|16|) and neglecting a 
q. (24), we finally 



sharp (S-peak due to the first term in 
obtain the power spectrum S(uj) of the electric field gen- 
erated by a resonantly-driven atomic system, 



S(u) = 
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of the atomic source. A closer look at Eqs. (19) and (20) 
shows that (E(x)) t is non-zero, even when the atomic 
system is in its ground state. When brought sufficiently 
close, the atomic dipole moment is expected to exert a 
force on a test charge. When this force is measured, the 
atomic state changes accordingly either into | Ai (£)) or 
into \X 2 (t)). We think that it would be interesting to 
observe this effect experimentally. 

IV. CONCLUSIONS 
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FIG. 3: The relative amplitudes of the Comparison of the un- 
normalised height of the maximum peak of S(u>) and Sm i(oj) 
for T = 10 8 Hz, ujo = 10 15 Hz and fl = 0.5 T. 



where the frequencies /i and 5 are as defined in Eqs. ([7]) 
and (12). This spectrum has some similarities with the 
Mollow spectrum S'moi(w) in Eq. (11) but also certain 



differences. 

Figs, [l] and [2] show the w-dependence of S (solid lines) 
and SmoI (dashed lines) for two different values of the 
laser Rabi frequency fi. Here the amplitudes of both 
spectra have been normalised to unity. When fl is rela- 
tively small there is only a single central peak. For suf- 
ficiently strong laser driving, both spectra exhibit three 
distinct peaks indicating both spectra contain similar in- 
formation about the atomic dynamics. There are, how- 
ever some differences. For example, the relative ampli- 
tudes of the sidebands of S(uj) are significantly smaller 
than the sidebands of Mollow's spectrum. 

Another difference between S(ui) and S'moi(w) is illus- 
trated in Fig. [3] This figure shows the dependence of 
the height of the central peaks with the laser Rabi fre- 
quency fl. The Mollow triplet is normalised such that 
JdojS'Moi(w) is a direct measure of the stationary state 
photon emission rate I ss of the laser-driven atomic sys- 
tem [TS1 [TH]. This means, its amplitude tends to zero, 
when fl tends to zero. In contrast to this, S(u) assumes 
its maximum when 51 — > 0. This means the spectrum 
S(u) is not a measure of the photon emission intensity 



This paper calculates the power spectrum S(u>) of the 
electric field (E(x)) t generated by a laser-driven atomic 
two-level system as would be observed at a relative posi- 
tion x. The detector should be placed a very small dis- 
tance away from the atomic source which is in principle 
feasible experimentally using currently available technol- 
ogy [6HTT]. The derived expression for S(uj) in Eq. (26) 



has some similarities with Mollow's resonance fluores- 
cence spectrum [T7HTO] but there are also several differ- 
ences. For sufficiently strong laser driving there is a cen- 
tral peak as well as two sidebands but their relative height 
is significantly reduced (cf. Figs, fl] and [2]). Moreover, the 
amplitude of this spectrum has a different dependence 
on the laser Rabi frequency fl (cf. Fig.|3|. It assumes its 
maximum when fl tends to zero. 

The main result of this paper is to identify the electric 
field generated by a laser-driven atomic two-level system. 
Due to the presence of a photon- absorbing environment, 
there is on average never more than one photon in the ra- 
diation field [T2TH5] , As a consequence, the only non-zero 
contribution to (E(x)) t comes from the atomic dipole 
moment r. This quantity can be calculated analytically 
without any approximations. The predicted spectrum 
S(lu) contains a similar amount of information about the 
atomic system dynamics as Mollow's spectrum but can 
be calculated in a more direct way and constitutes an 
interesting alternative property. We hope that this pa- 
per provides new insight into the dynamics of quantum 
optical systems with spontaneous photon emission. The 
present discussion can be extended relatively easily to 
more complex systems. 
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